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ABSTRACT 

The discovery that many trans-neptunian objects exist in pairs, or binaries, is proving 
invaluable for shedding light on the formation, evolution and structure of the outer 
Solar system. Based on recent systematic searches it has been estimated that up to 
10% of Kuiper-belt objects might be binaries. However, all examples discovered to- 
date are unusual, as compared to near-Earth and main-belt asteroid binaries, for their 
mass ratios of order unity and their large, eccentric orbits. In this article we propose a 
common dynamical origin for these compositional and orbital properties based on four- 
body simulations in the Hill approximation. Our calculations suggest that binaries are 
produced through the following chain of events: initially, long-lived quasi-bound bina- 
ries form by two bodies getting entangled in thin layers of dynamical chaos produced 
by solar tides within the Hill sphere. Next, energy transfer through gravitational scat- 
tering with a low-mass intruder nudges the binary into a nearby non-chaotic, stable 
zone of phase space. Finally, the binary hardens (loses energy) through a series of rela- 
tively gentle gravitational scattering encounters with further intruders. This produces 
binary orbits that are well fitted by Kepler ellipses. Dynamically, the overall process 
is strongly favored if the original quasi-bound binary contains comparable masses. We 
propose a simplified model of chaotic scattering to explain these results. Our findings 
suggest that the observed preference for roughly equal mass ratio binaries is probably 
a real effect; that is, it is not primarily due to an observational bias for widely sepa- 
rated, comparably bright objects. Nevertheless, we predict that a sizeable population 
of very unequal mass Kuiper-belt binaries is likely awaiting discovery. 

Key words: celestial mechanics - methods: N-body simulations - minor planets, 
asteroids - Kuiper Belt - binaries: general - scattering 



1 INTRODUCTION 

Binary systems occupy a special place in astronomy. First 
and foremost they allow the determination of binary partner 
masses which, in the case of stars, proved critical in the early 
development of the mass-luminosity relation iEddingtonl 
ll924l;lNoll2003h . In fact, the majority of all stars seem to be 
members of gravitational ly bound multiplets with most be- 
ing members of binaries (iHeggie fc Hutll2003r) . The dynam- 
ics and energetics of star clusters, in particular, are strongly 
influenced by the presence of binaries. This is because the 
largest single contributor to the total mechanical energy of 
the cluster may be the internal energy (binding energy) of 
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binaries llHillsil975l:[lnagakil984lJanesll99ll:lHeggie fc Hul 

2003). Not surprisingly the dynamics of binary systems and 



terest in stellar dynamics (lHeggieHl97Fl lHillslll975l 


1983allbl 


1990 


Mikkola & Valtonen 1992; Mardling 1995a b; 


Ouinlan 


1996 


Heggie & Hut 2003). 



Binaries are also relatively common in planetary 
physics; t hese include the Earth-Moo n and Pluto-Charon 
systems iChristy fc Har ringtonl Il978l : ICanup fc Asphaugl 
1200 it ICanupTT20M ~ l2bo5l) as "well as a substantial num- 
ber of asteroid binaries in the main belt. Asteroid bina- 
ries usually contain partners with mass ratios in the range 
m r ~ 10~ 3 — 10 -4 where m r = mijm\ and mi and m,2 are 
the masses of the primary and secondary binary partners 
respectively. Objects with such asymmetric m asses are most 
often thought of as "asteroids with satellites" iMerline et alJ 
2003) rather than as "binaries" although there exists no 
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Reference 

Path 1 <Weidenschillindl2002ft 



Path 3 iGoldreich et al.ll2002f> 



Path 4 iFunato et aI1l2004h 



Chaos-assisted capture (this article) 



Mechanism 

Physical collision and subsequent ac- 
cretion of two bodies inside Hill sphere 
of a third, larger Kuiper-belt object. 



Transitory binary formation inside 
Sun-binary Hill sphere. Capture and 
stabilization by dynamical friction. 



Transitory binary formation inside 
Sun-binary Hill sphere. Capture and 
stabilization by gravitational scatter- 
ing with a single large intruder - L 3 
mechanism. 

Initial production of asymmetric mass 
ratio "asteroid-like" binaries through 
physical collisions. Subsequent scat- 
tering encounters with large intruders 
then equalize mass ratio through ex- 
change reactions. 



Transitory binary formation 
through chao s - assis ted capture 
iAstakhov et alJ l200fl inside Sun- 
binary Hill sphere. Subsequent 
stabilization through a sequence of 
gravitational scattering encounters 
with small intuders each having ca. 
1 — 2% of total binary mass. 



Notes 

Depends on existence of approximately 
two orders of magnitude more massive 
bodies in the primordial Kuiper-belt 
than is currently accepted. At time of 
binary formation 99% of total mass 
of Kuiper-belt thoug ht to be in much 
smaller bodies jGgldreich et al ] boost 
iKenvon fc Bromlevll2004) . 

Predicts primordial binary formation 
rate ss 3- 10 — 6 per year per large body. 
Assumes larger sea of small (< 10 km) 
bodies than gra vitational instabilit y 
theories predict iFunato et alJ 2004). 
Predicts formation of multiplets. 

Predicts primordial binary formation 
rate s=s 3 ■ 10 — 7 per year per large 
body. No explicit mechanism provided 
for gradual reduction in post-capture 
binary orbit semimajor axes. 

Neglects solar tidal effects. Exclusively 
produces binaries with eccentricities 
e > 0.8. Effect of varying intruder 
mass not studied. Predicts future TNO 
binary discoveries will have roughly 
equal masses, large separations and 
high eccentricities. 

Includes solar tidal effects. Predicts ab 
initio the formation of (i) binaries with 
moderate eccentricities 0.2 < e < 0.8 
and (ii) a preference for roughly equal 
mass binaries. Also predicts that 
future higher resolution surveys of 
the Kuiper-belt will discover popu- 
lations of asymmetric mass binary 
TNOs. Mechanism suggests that 
binaries formed when the KB disk was 
dynamically cold, i.e., prior to excita- 
tion and depletion of the Kuiper-belt 
touncan fc Levisonll997UQomesl2003t 



Morbidclli . Emel'vanenko fc Levisorj 
20041) . 



Table 1. Proposed Kuiper-belt binary formation mechanisms. 



meaningful quantitative distinction between these defini- 
tions. The recent discovery of binary trans-neptunian ob- 
jects (TNOs) is particularly significant because these objects 
allow for the direct determination of TNO masses thereby 
opening up an important window into the origin, e volution 
and current composition of the Kuiper-belt (KB ) jStewart 
19971: lTothlll999l : iLuu fc Jewitl 120021: lNolJl2003t iGladman 
2005). It is quite remarkable that recent systematic searches 
iTruiillo fc Brownll2002l : lNol]||2003l: ISchaller fc BrowrJl2003D 
seem to suggest that up to 10% of K uiper-belt obj ects might 
exist as gravitationally bound pairs (Burns 2004). 

Kuiper-belt binaries (KBBs) are notable for their 
unusual compositional and orbital properties as com- 
pared to main-belt asteroids: the most salient are 
iNollI l200&t) (i) large mutual orbits, (ii) highly ec- 



centric mutual orbits and (Hi) m ass ratios of order 



unity (iMareotl 120021: 1 


Veillet et alJ 120021: iDurdal 


2002; 


Wcidcnschilline 2002: Goldrcich. Lithwick & Sari 2001 


llNoll 


2003: Schallcr & Brownl 


2003: lOsiD. Kern & Elliot 


20031: 


Altenhoff. Bertoldi & Mentenl 12004 IFunato et al.l 


12004 


Burns! 12004 iNoll et al. 


2004a b; Takahashi & Id 


|2004 


Nazzario & Hvda 120051). However, in each case certain 



caveats apply (iNollI 120031 ; INoll et"aH l2004allbT) . Firstly, by 



some measures KBB orbits are indeed large; e.g., for main- 
belt asteroids with satellites the semimajor axis of the mu- 
tual orbit (a) is typically about an order of magnitude larger 
than the radius of the pri mary asteroid (ta) whereas for 
KBBs a/r A ~ 50 - 500 iMerline et al.1 l2003t INoll et alJ 
2004b). Curiously, however, if the semimajor axis is com- 
pared to the radius of the binary's Hill sphere (Rh) then 
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main-belt asteroid binaries and KBBs have rather similar 
semimaj or axes, i. e ., expressed as a percentage a/Ru ~ 
2 - 5% (|NoUll2003l : INoll et al.ir2004bh . Secondly, although 
KBB orbits are noticeably eccentric as compared to the 
almost circular orbits of main-belt asteroid binaries, re- 
cent observations suggest that extremely large ec centricities 
(e > 0.8) may be rarer than was initially thought iNoli2003t 
INoll et al]l2004bft . Finally, while the mass ratios of known 
KBBs are in the range m r ~ 0.1 — 1 it is also possible that 
this result is at least partially due to an observation al bias 
for well-separated, equally bright objects feurnsl2004h . Nev- 
ertheless, the properties of KBBs are sufficiently striking to 
suggest that their formation mechanism differed consider- 
ably from that of main-belt or near-Earth asteroid binaries, 
which are thought to have been produced by physical colli- 
sions iCh auvmcau. Farinella fc Harrislll995l : IX u et al.lll995t 
iDurda et al.ll2005l) . 

Four different pathways (|W eidenschi llind 120021: 
iGoldreich et al.l l2002t IFunato et al.1 120041) - summarized 
in Table - have been proposed to explain the formation 
and properties of KBBs which we will refer to as Paths 
1-4: Path 1 involves physical collisions between two 
planete simals within the Hill sphere of a larger Kuiper-belt 
object iWeidenschillindl2002D . The two bodies then accrete 
and remain gravitationally bound as a single object to the 
larger body, thereby producing a binary. This mechanism 
assumes that binary objects are primordial; Sim ilar to 
Weidenschilling's mechanism \ Weidenschillinej I2002T) Paths 
2 and 3 also proceed from the temporary gravitational cap- 
ture of two objects within the three-body Hill sphere. Now, 
howe ver, the large third body is the Sun (IGoldreich et alJ 
2002). In both Paths 2 and 3 the initial stage of binary 
formation is the formation of a transitory binary which 
will eventually ionize unless stabilization occurs first. Two 
stabilization mechanisms have been proposed; in Path 2 
stabilization results from dynamical friction through inter- 
actions with a sea of smaller (< 10 km) Kuiper-belt bodies. 
In Path 3 the transitory binary is stabilized by energy loss 
through gravitational scattering with a large third body - 
the L 3 channel. Both of these mechanism s emphasize the 
importance of the Sun-b inary Hill sphere dSzebehelvlll967t 
iMurrav fc DermottllT999l) . 

A common feature of Paths 2 and 3 is that immedi- 
ately before and after capture the binary partners are fol- 
lowing an essentially three-body orbit. That is, solar pertur- 
bations or tides are import ant. Thus the orbits are periodic , 
quasi-periodic or chaotic jLichtenberg fc Liebermanl Il992l) 
and, in general, cannot be well described by orbit al elements 
of the Kepler problem (IMurrav fc Dermotlill999l) . However, 
actual KBB s follow orbits which can be well fitted by Ke- 
pler ellipses JVeillet et al]l2002t lNoJl2003t lOsip et aljl20o3 
INoll et al.ll2004allrJ) . i.e." solar tides are relatively unimpor- 
tant and the orbits are approximately two-body in nature. 
Therefore, a mechanism for gradually transforming captured 
- and, therefore, bound - quasiperiodic three-body orbits 
into (almost) periodic two-body orbits is required. We term 
this general process "Keplerization" - i.e., the gradual en- 
ergy loss of three-body binary orbits to produce orbits that 
are essentially Kepler ellipses. As three-body orbits lose en- 
ergy their semimajor axes undergo a steady reduction in 
size. In principle, Keplerization may occur directly during 
the capture process itself (e.g., as posited to happen in L 3 ) 



or it can happen more gradually after the initial capture 
event through, e.g ., continued dynamic al friction over ex- 
tended timescales <lGoldrei c^g^jlJ l2002l) . 

Path 4, prop osed bvlFun ato et all i2004l) . involves ex- 
change reactions (iHeggie fc HulJl2003l) in which the smaller 
member of an already bound "asteroid-like" binary is dis- 
placed by a larger body during a three-body encounter. This 
mechanism is essentially a v ari ant of binary star-intruder 
scattering <HillJll975L Il983allbl ll99(A iHut fc Bahca3ll99d : 
iHeerie fc Hutll2003D . In this mechanism solar tides are ig- 
nored and both the pre- and post-collision binary (if it sur- 
vives) follow a two-body Keplerian ellipse. 

While each of these four pathways is feasible, no con- 
sensus has yet emerged as to their relative importance or 
even whe ther a differe nt mechanism altogether might have 
operated iBurnsll2004h . This is because each mechanism ad- 
mits at least one potentially serious drawback. The first 
three paths are sensitive to assump tions about the siz e 
distributions of Kuiper-belt objects (IFunato et alJ 120041) . 
Path 1 depends on the existence of more large bodies in 
the K uiper-belt than seems to be consistent w ith observa- 
tions (IGoldreich et alJl2002l : IFunato et alJl2004h : The basis 
of Path 2 is dynamical friction which, to be effective, requires 
a much larger sea of small bodies th an is predicted by current 
theories of planetesimal formation iGoldreich fc Wardll973l: 
I Wetherhill fc Steward Il993f: iRafikovl 120031: iBernstein et alJ 
12004 IFunato et alJ 120041 : iKenvon fc Bromlevl 120041) . Fur^ 
ther, it is unclear how and why dynamical friction should 
select preferentially for roughly equal mass ratio binaries. 
This assumes, of course, that the apparent preference for 
mass ratios of order unity is not an observational artefact 

jBurnsll2004D . _ ^ 

The L 3 channel (IGoldreich et all2002l) . Path 3, depends 
on relatively rare close encounters between three large ob - 
jects and, based on the estimates of IGoldreich et all |2002), 
it is unclear if these occurred often enough to produce the 
current population of KBBs. In this regard it is interesting 
to note that recent calculations suggest that, if the Kuiper- 
belt lost its mass through collisional grinding, then an or- 
der of magnitude more binaries were likely present in the 
primordial Kuiper-bel t than has otherwise been thought 
iPetit fc Mousisl l20o4 - This illustrates that the number 
and size distributions of primord i al KB objects might 
be subject to revis ion ( Sterrjj2002£jBeraste^ne^L| 2004; 
[ Kenyon fc BromlevlEooi iPetit fc Mousdl2004l: lElliot etaU 
2005). A second drawback to Path 3 is that, as origi- 
nally proposed, it provides no explicit mechanism - be- 
yond possible post-encounter dynamical friction - for pro- 
ducing the approximately two-bo dy Keplerian binary or- 
bits which are actually obser ved (jyejnetetjrtJ |2002jjNoll| 
120031: I gchaller fc Brownl 120031 : lOsip et alJ 12003: INoll et*aTI 



2004a b). It is unlikely, as our simulations confirm, that a 
single collision between a quasi-bound binary (i.e., one that 
is following a three-body Hill trajectory) and a massive in- 
truder will result in a two-body Keplerian binary orbit. Fi- 
nally, Path 4 leads almost exclusively to binary eccentricities 
e ^ 0.8 a nd very large semim ajor axes which may approach 
Rh itself (IFunato et al.l20 04). However, moderate eccentric- 
ities, in the range 0.25 ^ e 0.82, and semimajo r axes that 
are only a few pe rcent of Rh seem to be the rule (lNollll2003l : 
lOsip et alJliool INoll et al.lEo04allrl . 

In the model described herein the first step is the 
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formation of long-lived quasi-bound binaries through the 
recentl y proposed mechanism of chaos-assisted capture 
fCAC. lAstakhov et alJ 1200.4 lAstakhov fc Farrellvl 12004 
Trim blel^AschwajTrien| |200j) in the Hill approximation 
JSzebehelvlll967l: iHenon fc Petitjll98^ : IMurrav fc Dermotd 
Il999l: iHeggie fc Hudl20olT In CAC particles become tem- 
porarily caught-up in thin chaotic layers within the Hill 
sphere of, in this case, the Sun-binary system. These chaotic 
layers, which separate regular from asymptotically hyper- 
bolic (scattering) orbits, are the result of the perturbation 
of the Sun on the motion of the two bodies making up the bi- 
nary. Because the orbits are chaotic any quasi-bound (tran- 
sitory) binary which is formed will eventually break apart. 
However, the lifetimes of these quasi-bound objects can be 
sufficiently long that, in the interim, stabiliz ation can occur. 
Since t he chaotic layers in the Hill problem iSimo fc Stuchil 
l2000l: lAstakhov et alJ 120031 : lAstakhov fc Farrellvl l2004fl 
lie adjacent to r egular K olmogor ov-Arnold-Mose r (KAM, 
IZaslavskvl Il985l : iLichtenberg fc Lieberman] I1992T) regions, 
then, in principle, even relatively weak perturbations can 
switch the binary into such KAM regions and thereby lead 
to permanent capture. 

Of course, the precise nature of the stabilization mech- 
anism is the crux of the matter. We propose that capture 
and subsequent hardening of thre e- body Hill orbits (i.e., 
Keplerization) jHillsl Il973 . Il983allbl Il99d: iQuinlanl Il996t 
IHeggie fc Hudl2003h proceed through chaotic gravitational 
scattering encounters with low-mass intruders which happen 
to transit the Hill sphere. Once a binary has been captured 
initially then further encounters with intr uders are proba- 
ble - about once every 3,000-10,000 years iWeidenschillinel 
2002) - and this can eventually produce binaries whose mu- 
tual orbits are Kepler ellipses. Empirically we find that the 
process is most efficient for equal mass binaries and rela- 
tively low-mass intruders having ~ 2% or less of the total 
binary mass. Incidentally, we argue that CAC is also the dy- 
na mical basis of the twin mechanisms proposed empirically 
bv lGoldreich et alJ l l2002h . 

To model the statistics of capture we compute cap- 
ture probabilities in four-body (Sun, binary, intruder) Monte 
Carlo scattering simulatio ns for various b inary mass ratios 
in the Hill approximation JScheeresi ri998). Our simulations 
reveal a clear preference for equal mass binaries to survive 
multiple subsequent encounters as they harden through in- 
truder scattering. 

The paper is organized as follows; Section [5] introduces 
the equations of motion for the three- and four-body prob- 
lem in the Hill approximation and provides a brief overview 
of how chaos can assist capture into stable three-bod y 
orbits jAstakhov et al.l l2003t lAstakhov fc Farrellvl I2004) . 
Mont e Carlo simula tions in the four-body Hill approxima- 
tion jScheeresUl99sl) are described in Sec. [3] A simplified 
(reduced dimensionality) model of chaotic low-mass intruder 
scattering is introduced in Sec. |^ to explain our results. To 
study the dynamics in this model we employ the Fast Lya- 
punov Indicator for several mass ratio and binary eccen- 
tricity combinations. A more detailed comparison with the 
other formation models summarized in Tabled is made in 
Sec.Qj] Conclusions are in Sec.[fj] 
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Figure 1. Typical level curves of the zero- velocity surface to- 
gether with the Lagrange po ints L\ and L2 for the thre e-body Hill 
problem in the planar limit iMurrav DermotJl99flft . The circle 
shaded grey shows the location of the Hill sphere. The Lagrange 
points demark the gateways into the capture zone. All units are 
scaled Hill units. 

2 EQUATIONS OF MOTION IN THE 

THREE-BODY AND FOUR-BODY HILL 
PROBLEMS 

In our model the initial stages of binary formation in- 
volve temporary capture within the Sun-binary Hill sphere 
through two particles getting entangled in chaotic regions 
of phase space, i.e., chaos-assiste d capture iAstakhov et alJ 
l2003t lAstakhov fc Farrellvl 1200^ . This process is well de- 
scribed in the three-body Hill approximation. 

2.1 Three-body Hill problem 

The three-body vector equ ations of motion in the Hill ap- 
proximation are given by (Szcbehclv Il967t IHenon fc Petid 
119861: IMurrav fc Dermotdll999t IScheeresHl998TT ~ 

p + Cl x [2p + n x p] = -p + 3a(ap) - (1) 

\P\ 

Here p — P2 — pi is the relative distance between the bi- 
nary members which have masses mi and mi and whose 
centre-of-mass is assumed to move along a circular orbit 
a = (1,0,0) around a larger body mo (the Sun). The co- 
ordinate system (x,y,z) is rotating with constant angular 
velocity CI = (0, 0, 1) which corresponds (in scaled units) to 
the mean motion of the binary centre of mass. In this coor- 
dinate system the centre-of-mass of the binary is located at 
the o rigin as is usual in the Hill problem jMurrav fc Dermotd 
1999). 

Units and typical values of param eters are defined a s 
follows with reference to 1998 WW31 iVeillet et alJl2002fl : 
mi ~ 2 ■ 10 21 g, ~ 10 21 p; the semimajor axis of the 
primary orbit at ~ 45AU and of the mutual binary orbit 
a ~ 22000 km. The Hill sphere (see Fig. occupies the re- 
gion between the saddle points L\ and L2 of the Hill problem 
which we term the "capture zone " iSzebehelvll967HScheeresl 
119981 IMurrav fc Dermotd Ir99sh : The Hill radius R H = 
a b ( m i+" 2 ) 1/3 ; in physical units: R H ~ 3 ■ 10 5 km ~ 25a. 
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In sc aled Hill units dSzebehelvl Il967l : iMurrav fc Dermotd 
I1999T) Rh = (1/3) 1//3 and the radii of the binary partners 



TA 
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2.2 Chaos-assisted capture in Hill's problem 

In order for a transitory (quasi-bound) binary to form two 
bodies must come inside their mutual Hill sphere defined 
with respect to the Sun. The Lagrange saddle points L\ and 
Z/2 in Fig. Q serve as gateways between the interior of the 
Hill sphere and heliocentric orbits. Figure^ shows a typical 
chaotic binary orbit obtained by integrating Hill's equations 
Q in two- dimensions. The orbit is trapped close to a peri- 
odic orbit for many periods before finally escaping from the 
Hill sphere (not shown) . The accompanying Poincare surface 
of section (SOS) in Fig. shows that the orbit is actually 
trapped in a chaotic layer separating a regular KAM island 
from a large region of hyperbolic scattering. Note that the 
orbits shown in Fig.|5]correspond to the mutual Hill orbit of 
the two binary partners and are indep endent of the particu- 
lar mass ratio of the bin ary partners jHenon fc Petidll986t 
IMurrav fc Dermotdll999l) . 

Exa mination of Poincare surfaces of section in the Hill 
problem dSimo fc Stuchill200d) - or, equivalently, the circu- 
lar restricted three-bod y problem (CRTBP) for small masses 
jAstakhov et al]l2003h - reveals that, at energies above the 
Lagrange saddle points L\ and L2 all phase space is di- 
vided into three parts, one of which regular KAM orbits 
inhabit, chaotic orbits another, the third consists of direct 
scattering orbits or, in other words, hyperbolic orbits. All 
these differ from each other in the behaviour of their or- 
bits. The chaotic orbits separate the regular from the hy- 
perbolic regions - see Fig.|2J). As energy is increased above 
the saddle points the chaotic regions consist of relatively 
thin layers of c haos which cling to the progressively shrink- 
ing KAM tori dSimo fc Stuchil 1200(3: lAstakhov et alJl2003t 
lAstakhov fc Farrellvll2003) . Because incoming orbits cannot 
penetrate the regular KAM tori in 2-dimensions (2D) and 
only exponentially slowly i n 3D llLichtenberg fc Liebermanl 
Il992l : lAstakhov et al.ll2003T) . orbits entering the Hill sphere 
from heliocentric orbits must either enter chaotic layers or 
scatter out of the Hill sphere essentially immediately. 

Binaries corresponding to chaotic orbits lying close to 
KAM tori can themselves follow almost regular orbits for 
very long times as shown in Fig. In the presence of a 
source of dissipation or of other perturbations these orbits 
can easily be switched into the nearby, but otherwise impen- 
etrable, KAM regions. 

It is remarkable that capture can occur if the binary de- 
creases or increases its energy during the encounter. This can 
happen because the topologies of stable KAM islands are, 
in general, a non-linear function of energy as is illustrated 
for four selected energies in Fig. That is, if a perturba- 
tion causes a chaotic orbit to increase or decrease its energy 
suddenly it may, in either event, find itself captured in a 
nearby regular part of phase space. Alternatively, such per- 
turbations, though small, may cause the quasi-bound binary 
to break up more quickly if the size of the chaotic layer is 
suddenly reduced in favour of scattering regions. 

Possible energy loss mechanisms which can sta- 
bilize transitory bi naries include physical collisions 
J Weidenschillind 12002ft . gravitational scattering with other 



bodies and dynamical friction dGoldreich et alJ I2002T) . 
Ruling out physical collisions and dynamical friction for 
reasons already described leaves gravitational scattering 
with "intruder" objects which enter the binary Hill sphere. 
The efficiency of this process depends critically on the 
intruder mass. For example, we find that the delicate 
chaotic binary orbits tend to be catastrophically desta- 
bilized by large intru ders as, e.g., in the L 3 mechanism 
jGoldreich et alJl2002lL Further, as Fig. |2]sliows, the orbits 
of transitory binaries can be a very large fraction of the 
Hill radius whereas actual binary orbits are typically only 
a few percent of Rh- Any KBB production mechanism 
must, therefore, explain how these large initial orbits are 
transmogrified into compact, Keplerian orbits. It is unlikely 
(as our simulations will show) that a single scattering 
encounter would propel the binary into such a final state 
directly. Rather, we propose that a series of gentle scattering 
encounters with low-mass intruders is necessary to capture 
and finally produce Keplerian orbits. Gentle stabilization 
by low mass intruder scattering is possible because only 
a weak perturbation is needed to drive long-lived chaotic 
orbits into adjacent regular regions of phase space. 



2.3 Four-body Hill problem 

The simulation of intruder scattering is facilitated by study- 
ing intruder scattering in t he four-body (Sun, binary, in- 
truder) Hill approximatio n iScheeresI Il998t iBrasserl I2OO2I : 
IScheeres fc Beilerosell2005ft . This approximation is appropri- 
ate given the mass ratios involved and the low to moderate 
eccentricities of typical KBB centre-of-mass orbits around 

the Sun. 

The generalization of the Hill problem l|Szebehelvll967l : 
iHenon fc Petidll986l: IMurrav fc Dermotdll999l) to the case 
where three small bodies, with a mutual centre-of-mass, R c , 
orbit a much larger fourth body (the Sun, mo = 1) on a near 
circular orbit (e.g., three Kuiper-belt obje cts within their 
mutual Hill sphere) is due to lScheerea il998T) . The total mass 
of the three bodies is defined by 



E 



rrij <S 1 



(2) 



where R c ~ a = (1,0,0) defines the motion of the three- 
body centre-of-mass along an almost circular orbit which 
defi nes the rotating frame. The vector equations of motion 
are JScheereslll99Sr) 

P 



p + ft x [2p + fi x p] = -p + 3a(a ■ p) - («i + a 2 ) 
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Here P3 is the coordinate of the third intruder body, 
m.3, and rrij = /zoy where 



3 



a, = 1. 



(5) 
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Figure 2. Non-linear dynamics of transitory binary formation in Hill's three-body problem. Panel (a) shows a quasi-bound (chaotic) 
binary centre-of-mass orbit (green) temporarily trapped close to a periodic orbit (red) which lies at the centre of a stable KAM island. 
This trajectory eventually escapes from the Sun- binary Hill sphere (large black circle). Panel (b) is a composite showing the colour coded 
SOS [x — y hyperplane defined such that the momentum p x = and the velocity -j| > 0) for the trajectories in panel (a) and also several 
nearby chaotic (grey) and quasi-periodic (blue) orbits. Panel (c) shows the non-linear evolution of the size and shape of the KAM island 
in (b) with energy at four selected energy values. Distances and energies are in scaled Hill units. 



Equations (J^J and Q contain the coordinates pi,p2 
explicitly; this is purely for compactness of notation and, in 
practice, they are eliminated using the centre-of-mass rela- 
tion 

3 

^2<XjPj=0. (6) 

3=1 

When 7Ti3 = equation @ re duces to the three- 
body Hill problem Q f Szebehervll967t iMurrav fc DermottJ 
1999) and becomes uncoupled from equation Q. Equa- 
tion © i s th en a var i ant o f the restricted Hill four- 
body problem JScheeresI Il99gl) and depends explicitly on 
the solutio n p(t) which may be periodic , quasi-periodic, 
or ch aotic llLichtenberg fc Liebermanlll992l lAstakhov et alJ 
2003). These solutions are actually particular trajectories 
for the three-body Hill problem. 

To model the statistics of capture immediately after 
transitory binary formation we computed capture proba- 
bilities in four-body Monte Carlo scattering simulations as 
described in the next section. 



3 MONTE CARLO SCATTERING 
SIMULATIONS 

We are primarily interested in binary-intruder scattering, 
with intruders approaching from infinity. Therefore we as- 
sume that the transitory binary is initially following a 
chaotic, but long-lived three-body Hill orbit. 

Because our Monte Carlo scattering simulations involve 
chaotic binary orbits they differ in essential respects from 
previous simulations; to our knowledge all previous simu- 
lations have dealt with intruder scattering from bound bi- 
naries whic h follow circular or elliptical two-body Keple- 
rian orbits lHiUdll975L Il983allbl Il990t [Hut fc Bahcalllll99ri 
iHeggie fc Hutj2003l : lFunato et alJl2004T) . Therefore it is nec- 
essary to explain our procedure in detail. 

In any energy range for which chaotic layers exist it is 



possible to find an infinite number of initial conditions cor- 
resp onding to different quasi-bound, three-body binary or- 
bits llLichtenberg fc Liebermanlll992l) . This is illustrated in 
Flg-Ht; it is apparent that a swath of chaotic orbits over the 
energy range shown will not only have a range of lifetimes 
but their dynamics will depend sensitively on their partic- 
ular initial conditions and energy. Since we are interested 
in the probability of capture from chaotic zones into stable 
KAM zones it is, therefore, necessary to study an ensem- 
ble of such orbits rather than a single example. That is, the 
capture probability of interest is not simply the probability 
of any given orbit being captured; it must also measure the 
density of "capturable" orbits inhabiting the chaotic zones. 

3.1 Methods 

The procedure we have developed is as follows: 

(i) Initially a cohort of 170 different quasi-bound chaotic 
binary orbit initial conditions was harve sted by int egrat- 
ing eq. <|T} usi ng a Bulirsch-Stoer method jPress et alJll999l : 
lAarsethll2003ri . Initial conditions were chosen randomly in- 
side the three-body (Sun-binary) Hill sphere for randomly 
chosen negative energies lying above the (three-body Hill) 
Lagrange saddle points and in the energy range —2.15 < 
E < —2.0. If an orbit remained inside the Hill sphere for a 
time Ti < Th < T2 then its initial conditions were stored. 
We refer to Th as the orbit's "natural" Hill lifetime. In scaled 
units we chose Ti = 15 (~ 600 years) and Tb = 50. While 
much longer (and shorter) living orbits are possible selected 
integrations for such orbits produce comparable results. This 
procedure guarantees that all orbits in the cohort are chaotic 
since, absent perturbations, they all eventually exit the Hill 
sphere. 

(ii) Equations and were integrated for various val- 
ues of m r and 7713 using the stored initial conditions for the 
binary, but rescaled to the actual values of mi and 1712 used. 
Thus, the same cohort of three-body quasi-bound Hill orbits 
was used for all mass ratios. 
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(iii) Intruder particles were fired at each binary orbit with 
the phase of the binary chosen randomly. For each binary 
orbit 5000 intruders were sent in and, each time, it was 
recorded whether the binary was stabilized or not. For each 
binary orbit this produces a capture probability, i.e., the 
fraction of the 5000 intruders that lead to stabilization in a 
single encounter with a binary. 

(iv) Intruders were started isotropically with uniformly 
chosen random positions on a sphere of radius Rh, with 
uniform rando m velocities \v \ ^ 5 v h where vh is the 
Hill velocity jGoldreich et alJ |2002t iFunato et all |2004 
lOoldreich. Lithwick fc Sa.rill20ol~ 

v h ~£IRh, (7) 

with being the orbital frequency of the third body around 
the binary at distance Rh- 

(v) In order to choose the phase of the binary randomly 
intruder integrations were started at randomly chosen times 
t G (0, Th) along the binary orbit. In practice this has the 
effect of sending the intruder towards the binary at different 
relative configurations of the binary partners in phase space. 

(vi) The integrations proceeded until either the binary 
broke up or it survived for a time t — 10 Th at which 
point it was considered captured. After the intruder escaped 
from a sphere of radius 3Rh the integrations for the binary 
continued assuming = 0, i.e., the usual three-body Hill 
equations Q were used. 

(vii) At the end of the integration each captured binary 
and its intruder were back integrated in time to ensure that 
the binary (in the four-body integration) broke up in the 
past. This helps protect against accidentally starting the 
binary-intruder trinary in already bound regions of phase 
space. 

(viii) Integrations were stopped if particles came within 
a distance ta = 10 -4 of each other. 

(ix) These simulations were repeated ten times so as to 
be able to compute variances (error bars) for the capture 
probabilities. Each time a different random number seed was 
used to guarantee that the simulations were independent. 
Thus, e.g., Orbit 1 may end up with capture probability 
0.4 ± 0.02, Orbit 2, 0.23 ± 0.018, and so on. 

(x) Orbits in the cohort were then binned according 
to their assigned capture probabilities. This produced his- 
tograms showing the number-density of orbits in the cohort 
having capture probabilities in the various specified inter- 
vals. Mainly the error bars were used (a) to establish that 
the calculations had converged and (b) to decide on a rea- 
sonable granularity for the probability bins. 

(xi) The actual number of binary orbits included in the 
cohort was chosen so as to provide acceptable convergence 
given available computational resources. Convergence was 
judged acceptable when variances in capture probabilities 
were on the order of a few percent or less. 



3.2 Single intruder scattering 

In these simulations, each binary orbit in the cohort ends 
up with a capture probability and an error bar (of the order 
of a few percent) attached to it. However, for different mass 
ratios the same binary orbit in the cohort will, in general, 
have different capture probabilities and error bars. Figure 
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Figure 3. Fraction of binary orbits (number-density distribu- 
tions) binned according to their individual capture probabilities 
as described in the text with 03 = 0.025; Mass ratios are (a) 
a\ = 02 and (b) ol\ = 14c*2 (b). Smooth curves are best fits 
to normal distributions. The same set of 170 chaotic Hill orbits 
in three-body centre-of-mass coordinates, appropriately rescaled, 
was used in both cases. The binary was considered captured only 
if the original binary partners remained bound to each other, i.e., 
exchange reactions were ignored. In 98.5% of cases the average 
capture probability for any given Hill orbit was higher for equal 
mass binaries. 

presents histograms showing the fraction of captured orbits 
binned according to their individual capture probabililities 
for two typical mass ratios. We tested that these histograms 
are robust to the actual number of binary orbits included in 
our cohort - e.g., if 85 rather than 170 orbits are used the 
shapes of the histograms are essentially identical to those in 
Fig. |3] The width of the bins is roughly commensurate with 
the average width of the error bars. 

The distributions in Fig.[3]clearly indicate that the frac- 
tion of binary orbits peaks at a significantly higher capture 
probability for equal mass ratio binaries than for unequal 
mass ratio binaries. This translates into a higher probability 
for capture in the equal mass case. Recall that the actual bi- 
nary orbits used in both cases are identical. Since the cohort 
of binary orbits used is the same, and the intruder masses 
and initial conditions are also the same, then the difference 
in the distributions is due entirely to the different binary 
mass ratios employed. In fact the only parameter that was 
varied between the simulations shown in Fig.|3]was the mass 
ratio. 

We also compared individual capture probabilities or- 
bit for orbit and found that in 98.5% of cases the aver- 
age capture probability for any given Hill orbit was higher 
in the equal mass simulation. Thus, at any fixed Hill ra- 
dius (i.e., the same total binary mass) a clear preference 
for the capture of same-sized binary partners is appar- 
ent. For comparison, the mass ratios m r of the four best 
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chara cterized KBBs are ~ 0.57 ( 199 8 WW™ . IVeillet et alJ 
12002ft . ~ 0.56 ((66652) 19 99 RZ 2S 3, [Noll et al.ll2004al) . ~ 
0.55 ((58534) 1997 C O™. iNoll et all l2004rjT and ~ 0.34 
((88611) 2001 QT 297 , lOsip et alj|2003h . These are signifi- 
cantly larger than the mass ratios o f main-belt binaries for 
which, as noted, m r ~ 10" 3 - 10" 4 llMerline et al.ll2003D . 

The calculations for Fig. [3] took ~ 1 month using all 
nodes on a 32-node Beowulf cluster. Interestingly the mass 
effect manifested itself quite directly in that the simula- 
tions for Fig. [3p finished about a week earlier than those 
for Fig. This is because of the additional tests and inte- 
grations that needed to be performed to verify that an orbit 
had been captured permanently. 

3.3 Multiple intruder scattering 

The probabilities in Fig. [3] are for permanently captur- 
ing already-formed quasi-bound binaries through a sin- 
gle intruder scattering event. However, the overall prob- 
ability of binary formation depends too on the probabil- 
ity of forming quasi-bound binaries in the first place. Be- 
cau se Hill's equation s ffl are parameter free (after rescal- 
ing. iMurrav fc Dermotllig gg 1 *) once two objects come within 
their mutual Hill sphere the probability of chaos-assisted 
quasi-binary formation becomes independent of their rela- 
tive masses. Nevertheless, the entry rate of bodies into each 
other's Hill sphere - and thus the overall capture proba- 
bility - will depend on the actual masses and velocities of 
the objects involved. Uncertainties in the size and velocity 
distributions of contemporary and primordial Kuiper-belt 
objects means, however, that estimating Hill sphere entry 
rates is an open problem but one whose resolution should 
be aided by t he discovery and characterization of further 
KBB objects JPetit fc Mousd 12004 iBernstein et aLlliool 
iKenvon fc Bromlevl 120041: INoll et al.ll2004b() . Despite these 
uncertainties the mass effect is expected to persist because 
it is progressively enhanced by later scattering events. 

After its initial capture through a single scattering en- 
counter with an intruder a binary is, typically, following a 
very large three-body orbit with a semimajor axis similar in 
size to that of the periodic orbit illustrated in Fig. That 
is, though now bound, the binary partners are still strongly 
influenced by solar tides. In this subsection we demonstrate 
that essentially Keplerian two-body orbits can be produced 
through a series of subsequent scattering encounters with 
further small intruders. 

That this is possible is demonstrated by the simulations 
underlying Fig. 2] In these simulations, for each mass ratio, a 
cohort of 1200 randomly chosen Hill binary orbits (in the en- 
ergy range —2.15 < E < —0.15) was generated as described 
in the previous subsection. The scatter plot shows the final 
orbital elements of only those binaries out of the 1200 which 
survived N enc = 200 successive encounters with intruders. 
Intruder initial conditions were generated as described pre- 
viously. In this set of simulations an encounter was consid- 
ered non-disruptive if the binary survived for t > 50 scaled 
Hill time units. The first intruder was sent in at a random 
point along the binary orbit during the orbit's Hill lifetime 
t £ (0, Th) (the binary orbit integration was itself started at 
t = 0) . If the binary survived to t = 50 the next intruder was 
sent in. If the binary survived that event then the next in- 
truder was sent in, and so on. Thus, intruders were spaced 50 
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Figure 4. Distributions of orbital elements, eccentricity, e, and 
semimajor axis, a (scaled by the Hill radius Rh) after N Enc = 200 
intruder scattering events per binary The figure combines results 
for twenty mass ratios: m r = 22. = 1 (blue); m r = 0.05 (red) and 
18 intermediate values 0.1 Sj m r ^ 0.95 (green dots). Yellow dia- 
monds show the locations of the four KBBs with known orbital el- 
ements (boxes are la observational er ror bars) 1 



/eillet et al.l2l 

lOsip et al.ll2003t INoll et alJl2004alrl . For each mass ratio, 1200 
different initial binaries were integrated to generate the figure. 
To accumulate statistics this simulation was repeated 60 times for 
the two extreme mass ratios. In each case 5000 (rather than 1200) 
initial binaries were used per mass ratio. After N enc = 200 en- 
counters binary survival probabilities were 0.103±0.007 (m r = 1) 
and 0.019±0.002 (m r = 0.05). The survival probability is defined 
as the fraction of the orginal 5000 binaries which remains bound 
after N enc = 200. 



time units apart. This was repeated until either the binary 
broke up or it survived N enc such successive encounters. 

In order to generate statistics it proved desirable to per- 
form a much larger simulation than was actually used to pro- 
duce Fig. 0] The capture probabilities and error bars given 
in the caption to Fig. ^Jwere obtained by repeating the sim- 
ulations as described in 60 independent runs, in each case 
using an initial cohort of 5000 (instead of 1200) randomly 
chosen binary orbits. 

Naturally, each of the binary-intruder encounters may 
individually be stabilizing or destabilizing. However, our 
simulations indicate that comparable mass binaries, es- 
pecially, tend to survive multiple encounters as they 
harden through intruder scattering, thereby becoming more 
Keplerian; this is similar to the situation in low-mass 
intrudcr/binary-star scattering in which (depending o n in- 
truder veloc it y) orbits tend, on averag e, to harden jHillsl 
ll975lll983allblll990l : lHeggie fc Hutl2003l) . After N enc = 200 
intruder encounters Fig. 2] demonstrates that multiple scat- 
tering can produce KBBs with orbital elements comparable 
to those actually observed. Once a binary has been captured 
then further encounters with intruders are likely. We assume 
~ 100 km-sized binary partners and intruders having masses 
1 — 2% of the total binary mass. Based on Weidenschilling's 
model we find that binary-intruder encounters will occu r 
about once every 3,000-10,000 years <IWeidenschillindl2002Tl : 
Fig-Bldemonstrates that such encounters can eventually pro- 
duce Keplerian orbits. Although the final orbits in Fig.0|arc 
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still strictly solutions to Hill's problem, their energy has been 
reduced sufficiently that, for all intents and purposes, they 
follow two-body Kepler ellipses, i.e., th e orbits are well de- 
scrib ed by Keplerian orbital elements jMurrav fc Dermotd 
Il999ft . 

Experimental simulations for multiple intruder encoun- 
ters reveal that the main effect of increasing the number of 
encounters (N enc ) is to reduce the final value of the semi- 
major axis whereas the final eccentricity distributions were 
quite robust to N enc . Thus we chose the single parameter 
N enc = 200 heuristically so that orbits ended up roughly in 
the observed semimajor axis range of 1998 WW31. There- 
fore, Fig. ^represents a single parameter fit to the semimajor 
axis of one of the most well characterized KBBs. However, 
no special attempt was made to refine the fit. 

Kuiper-belt binary semimajor axes - expressed as a 
fraction of the binary Hill radius - are generally compa- 
rable to those of asteroid binaries and Fig. [I] captures this 
finding. Figure 2] also reproduces the moderate eccentrici- 
ties which seem to be a feature of actual KBBs jNoll et alJ 
l2004bt) . This finding contrasts with the extremely large ec- 
centricit ies (e > 0.8) produ ced exclusively in exchange re- 
actions jFunato et al ■1120041) . Finally, Fig. 31 demonstrates 
that equal mass binaries have significantly higher survival 
probabilities than asymmetric mass binaries after multiple 
intruder encounters. 



4 REDUCED MODEL OF INTRUDER 
SCATTERING 

The simulations described so far indicate that roughly equal 
mass binaries have a statistically significantly higher cap- 
ture probability than binaries with small mass ratios. How- 
ever, the origin of this effect is not directly apparent from 
the simulations. In this section we develop a reduced model 
of intruder scattering which provides an explanation of the 
mass effect discovered in the simulations. The origin of this 
effect is related to chaotic scattering of the intruder by the 
binary. 



4.1 Chaotic scattering and dwell times 

Chaotic scattering involves a complex and sensitive depen- 
dence of some "output variable" (e.g., scattering angle, 
interplay- or dwell-time, etc.) on an "input variable" (e.g., 
impact parameter). For example, scatter ing in the CR TBP 
has been demonst r ated t o be chaotic jBoyd fc McM illan 
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iMacau fc Caldasl 120021: iNaelerl 

plotted scattering functions (out- 
put energy and trajectory length) as a function of impact 
parameter and found clear evidence of chaotic scatte ring. 
Similar results were found bv lBovd fc McMillanl <1993ft who 
demonstrated chaotic scattering in binary star - intruder 
scattering. However, the high dimensionality of the four- 
body problem, even in the Hill approximation, makes it dif- 
ficult to study a single output variable as a function of a 
single input variable. Nevertheless, we performed a series of 
exploratory simulations in which the dwell-time of the in- 
truder inside the Hill sphere was calculated as a function of 
binary mass ratio. The dwell-time is here defined to be the 
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Figure 5. Dwell-time distributions for intruder scattering for 
(a) m r = 1 and (b) m r = 0.035. For each of the two mass ratios 
shown 5000 low-mass intruder scattering events were simulated 
as described in Sec. [3] The histograms show densities of intruders 
in various dwell-time ranges. Some dwell-times t > 10 occur but 
these have been omitted from the figure for clarity. For the sake of 
comparison, in these simulations, a single chaotic binary orbit was 
used but the results are similar if different chaotic binary orbits 
are used. An identical series of intruders was used in both cases. 
The Hill lifetime of the binary orbit used is Tjj 17 in scaled 
units. Intruders are binned according to the time they remain 
within the Hill sphere — their dwell-time. Average dwell-times are 
(a) 2.1 and (b) 2.7 scaled time units. 



total time the intruder remains inside the Hill sphere in our 
Monte Carlo single-encounter scattering simulations. 

Fig. |^| presents typical histograms of intruder dwell 
times for two mass ratios. The dwell-time distributions for 
the two cases differ in that for unequal masses the dis- 
tribution has a longer tail corresponding to longer dwell- 
times. That is, the intruder has a gre ater chance of b ecoming 
caught up in a long-lived resonance (iHeggie fc Hudl2003ft if 
the binary has very asymmetric mass partners. The average 
dwell times shown are quite typical and are fairly consistent 
between different binary orbits. The unequal mass case typi- 
cally has an average dwell-time w 20 — 30% longer than that 
for equal masses depending on the actual mass ratio used 
and the particular binary orbit selected. 

In an attempt to understand this result and how (or if) 
it relates to the differential stabilization of equal and un- 
equal mass binaries we have developed a simplified model 
of intruder scattering. The main reason why we resort to 
studying a simplified problem is the difficulty in visualizing 
the dynamics and the structure of the 12- dimensional phase 
space in the full four-body Hill problem. We stress, however, 
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that, in the following, we are only interpreting the dynamics 
exhibited in the full problem using reduced dimensionality 
dynamics. None of the results obtained in the previous sec- 
tion rely on this analysis. 

Empirically it seems likely that the origin of the mass ef- 
fect in Figs.03and|I| can be traced to the observed difference 
in in truder dwell- or interplay-times ijShebalin fc Tippensl 
1996). This is because the only consistent difference we were 
able to find between simulations which were otherwise iden- 
tical, except for having different mass ratios, was in intruder 
dwell-time distributions. A similar effect has also been ob- 
served by Hills in s im ulatio ns of star stellar-binary scattering 
(|HmsJll97fll Il983allbl fl99flh . However, Hills took the binary 
orbits to be bound, two-body Keplerian orbits rather than 
chaotic three-body orbits; i.e., Hills' work does not contain 
the equivalent of solar tides which, in our model, are essen- 
tial for the formation of quasi-bound binaries in the first 
place. 

The model we build has its foundations in the following 
observations; 

(i) A separation in timescales exists between fast intruder 
scattering and the timescale of the mutual binary orbit. That 
is, chaos in the binary orbit caused by solar tides generally 
develops on a slower timescale than does the intruder scat- 
tering event. 

(ii) The difference in dwell-times between equal and un- 
equal mass binaries is greatest for small intruders. 

(iii) If simulations are done using zero-mass intruders 
then, clearly, no binary stabilization can occur. Neverthe- 
less, the intruder dwell-time effect persists. In this limit the 
problem reduces to the so-called restricted Hill four-body 
problem (RH4BP) (IScheereslll99il) . 

(iv) Experimental simulations reveal that the dwell-time 
effect in the RH4BP exists even if the binary follows an el- 
liptical Keplerian orbit. These findings are based on exten- 
sive simulations in the four-body Hill problem and in several 
variants of the RH4BP in which different solutions for p(t) 
were used - see the discussion following equation ©. 

(v) As the binary hardens then solar tides become rela- 
tively less important and the problem red uces to the ellip tic 
restricted three-body problem (ERTBP) (IScheereslll9"98n . 



4.2 The Elliptic Restricted Three-body Problem 

Combining the observations in the previous subsection leads 
to the following set of assumptions: neglect solar tides, as- 
sume elliptical binary orbits and set the intruder mass to 
zero. After making these approximations the four-body Hill 
prob lem reduces to the elliptic restricted three-body prob- 
lem llSc heeresll l998t) for which the Hamiltonian is given by 
JSzebeherylll967t iLlibre fc Pifioll ll99oL iMikkola fc Valtonenl 
ll992l:lAstakhov fc Farrelryll2004D ; 



1 2 



2 . 2\ 
Pz + Z ) 



+ ■ 



-^/(l + x) 2 + y 2 + z 2 ^Jx 2 +y 2 + z 



l (x 2 +y 2 + z 2 ) + (1 _^ )x+ l (1 _ M , 



= + 



Here the semimajor axis of the mutual binary orbit a has 
been scaled to unity; e is the eccentricity of the binary or- 
bit and fi' = m^j(m\ + ni2) or m r = ///(l — I 1 )• The true 
anomaly of the binary mutual orbit /, i.e., the angular posi- 
tion measured f rom the pericent re, is related to the physical 
time t through JSzebehelvlll967f) 
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The CRTBP Hamiltonian jAstakhov et alj [2003) 
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is recovered when e = 0. The Hamiltonians ||HJ and 1101 are 
defined in the rotating frame with the origin transformed to 
one of the primaries. 

Note that, in the following, the eccentricity being re- 
ferred to in the context of the ERTBP is the eccentricity of 
the mutual orbit of the binary partners and not the helio- 
centric eccentricity of the binary centre-of-mass. The binary 
centre-of-mass is still assumed to be following a circular he- 
liocentric orbit. 

It is important to keep in mind that, in the ERTBP, 
the two primaries are assumed to be following an elliptical 
(Keplerian) mutual orbit and, therefore, these bodies cannot 
ionize. In our simplified model the ERTBP primaries are 
actually the binary partners which are, in reality, following 
an unbound chaotic orbit. However, if the intruder dwell- 
time is sufficiently small then, during the intruder scattering 
event itself, neither the intruder nor the binary will sense the 
longer timescale chaos caused by solar perturbations on the 
binary orbit. In a sense this picture is similar to the Born- 
Oppenheimer approximation for molecules in which the slow 
nuclei appear almost frozen from the point o f view of the 
much faster electrons jKarplus fc Porterlll97Gft . It should be 
emphasized that, even if the intruder does not sense that the 
binary orbit is chaotic, the perturbation the intruder induces 
on the binary can, nevertheless, influence the binary orbit's 
long time development. That is the intruder can stabilize, 
destabilize, or leave essentially unaffected the evolution of 
the binary mutual orbit. 

The ERTBP will, therefore, only be a reasonable model 
if the timescale of the binary partner relative motion is much 
slower than that of binary-intruder scattering encounters. 
If this is the case then the four-body Hill problem can be 
approximately decomposed into two simpler three-body sys- 
tems; Sun-binary and binary-intruder. In practice, we find 
that intruder dwell times - the time the intruder spends in- 
side the Hill sphere - in the full four-body simulations are 
normally no more than a few percent of the Hill lifetime of 
the orbit (Th, defined earlier) and so we expect that study- 
ing the reduced (ERTBP) dynamics will be useful as a proxy 
for the full problem. 

Stated somewhat differently, we make the assumption 
that the relative approach velocity of the massive binary 
partners is quite small whereas hardening requires that the 
intruder approach in a hyperbolic flyby at a somewhat 
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larger relative velocity. While this might seem contradic- 
tory, equipartition of energy in a particle disk implies that 
the largest bodies will have the most circular coplanar or- 
bits, while the smaller bodies - here the intruders - will 
be more eccentric and inclined. Thus the relative veloc- 
ity of a pproach should indeed be larger for the small in- 
truders jGoldreich et all 12001 IWetherhill fc Stewardll993t 
Stewart fc Idall200Ct iLewis fc Steward 120021: lOhtsuki etld] 
20021: iGoldreich et aljl2004l . 

Because the binary is actually following an unbound 
chaotic orbit it is easily disrupted. Even small intruders can, 
in principle, lead to early disruption of the binary (as com- 
pared to Th, the binary lifetime in the absence of intrud- 
ers). If, on the one hand, the intruder has a short dwell-time 
then it is in and out of the Hill sphere quickly and there 
is an opportunity for sudden energy transfer with minimal 
disruption of the binary. Granted, this energy transfer may 
destabilize the binary but, in the intruder velocity regime 
of interest, we find that, on average, this is not the c ase - 
again this re su lt is q uite similar to the findings of Hills (Hills 
Il975l.ll983allbl 1990). Note especially that ionization of the 
binary is possible - though not necessarily equally likely - 
for small as well as large intruders. This situation should be 
contrasted with the case of already bound (Keplerian) bina- 
ries when, in general, a comparably massive intruder may 
be needed to produc e complete or "democratic" ionization 
jHeggie fc Hudl2003D . 

On the other hand if the intruder gets trapped in a rel- 
atively long-living resonance with the primaries then there 
exists a greater opportunity for disruption since the pertur- 
bation to the binary orbit is being applied for longer. That 
is, if the intruder gets caught up in a sticky chaotic layer 
its stay within the Hill sphere will be extended (in a reso- 
nance) as compared to if it had entered an asymptotically 
hyperbolic regime. While the lifetime of the resonance may 
still be much shorter than Th it will be significantly longer 
than the dwell time of an intruder undergoing hyperbolic 
scattering. Such resonances have t hree main decay channels 
jHilldll983bl iHeggie fc Hudl2003fl : (i) complete ionization, 
(ii) expulsion of the intruder and (iii) expulsion of one of 
the original binary partners (exchange). While channel (ii) 
is the most likely of the three for small intruders, channels 
(i) and (iii) can also occur because the binary orbit is not 
bound. Thus, as compared to sudden scatteri ng, resonances 
(in th e sense of long-lived trinary complexes iHeggie fc Hud 
2003)) tend, on average, to reduce the capture probability. 

In view of the foregoing we suggest that the relative 
sizes of the chaotic regions as compared to the hyperbolic re- 
gions will correlate directly with the intruder's dwell-time. 
We argue that the relative sizes of these regions should, 
therefore, also correlate to capture probabilities. Chaos de- 
lays the intruder within the Hill sphere and so amplifies the 
effect of the intruder on the binary orbit. Repeated inter- 
actions - and energy transfer - between the binary and the 
intruder tend to destabilize the fragile binary orbit. Impor- 
tantly, we are here comparing relative rather than absolute 
propensities for stabilization. Not all fast encounters are sta- 
bilizing nor do all resonances cause the binary to decay. 

Thus, what requires examination is the relative size 
of chaotic versus hyperbolic regimes for various mass ra- 
tio/eccentricity combinations. The value of the ERTBP 
limit of the four-body Hill problem is that it allows one 



to visualize phase space directly; in this case using the 
Fast Lyapunov Indicator (FLI) as will now be describe d 
jFroeschle. Guzzo fc Legal200d : lAstakhov fc Farrellvl2004h . 

4.3 Fast Lyapunov Indicator 

The FLI is useful for mapping chaotic and regular regions of 



phase space when Poincare SOS cannot readily be computed 


jFroeschle et all 1200(1 


Pilat-Lohinger. Funk & Dvorak! 


20031 lAstakhov & Farrellv 


2004). Given an n-dimensional 



continuous-time dynamical system, 



dx/dt = F(x,i),x = (X!,X 2 , ...,£n), (11) 

the F ast Lyapunov Indicator is defined as iFroeschle et all 
2000) 

FL/(x(0),v(0),t)=ln|v(i)|, (12) 

where v(t) is a solution of the system of variational equa- 
tions iTancredi. Sanchez fc RoidboOll) 

dt = {dx-) V - (13) 

Regularization IStiefel fc Scheifeldll97ll:lAarsethll2003T) 
was used to integrate equations dill and (1131) . 

4-3.1 Choice of initial conditions and energies in FLI 
maps 

To analyze the structure of phase space in our simplified 
system, i.e., in the ERTBP, we computed FLI maps in the 
planar (2D) circular and elliptical cases. In the circular limit 
a direct comparison with surfaces of section (in the Hill limit, 
fi' < 1. lAstakhov fc Farrellv! I2004T) can be made. We have 
previously demonstrated that in the planar CRTBP the FLI 
reproduces correctly all rel evant features visible in the SOS 
llAstakhov fc Farrellvl2004T) . In particular, the chaotic layer, 
as it evolves with increasing energy, can be easily identi- 
fied by high values of the FLI. Note also that FLI measure- 
ments make sense even for relatively short time integrations, 
and, even in 2D, are much more efficient than is the con- 
struction of SOS flFroeschle et^lj 2000l ; |Pilat-LoIiinger et alJ 
120031 lAstakhov fc Farrellvll2004h .~ 

In the planar limit of the ERTBP initial conditions were 
generated randomly within the Hill radius on the surface 
of section (see Fig. [5] caption) taking, initially, / = 7r/2. 
This guarantees t hat the E RTBP initial conditions reduce to 
those of CRTBP (lAstakhov fc Farrellvll2004ll . For a given (i 
and e all ERTBP initial condi tions are, therefore, generat ed 
with identical initial energies (lAstakhov fc Farrellvll2004T) . 

A further issue in generating the FLI maps in Fig. HJ 
concerns the choice of initial energy. Because we are com- 
paring different masses and ellipticities it is not possible to 
work at exactly the same energy in each case. This is be- 
cause in the ERTBP the energies of the Lagrange points, 
for example, depend on ellipticity, mass and the instanta- 
neous value of t he true anomaly, /. That is, energy (or the 
Jacobi constant, iMurrav fc Dermotdll999h is not conserved 
in the ERTBP. Thus, in order to make a fair comparison of 
the dynamics between these different cases we need to con- 
struct a criterion of comparability; in practice, we selected 
initial values of the Jacobi constant such that the sizes of 
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Figure 6. Fast Lyapunov indicator maps for the planar ERTBP and four eccentricity/mass ratio combinations. The origin is shifted so 
that the binary partners are located at (x,y) = (±0.5,0) in scaled units and here units are rescaled so that Rjj = 1. Colour scale runs 
from blue to red to white with increasing FLI. Initial energies are as defined in the text. Blank regions correspond to directly scattering 
or hyperbolic tra je ctories which survive for les s than ~ 20 binary periods. Areas that resemble shotgun-like patterns of points correspond 
to chaotic regions I Asta khov &: Farrellvll^OO^ (see the text). At the higher eccentricities scattering (hyperbolic) regions are noticeably 
larger for (c) (equal masses) than for (d) (unequal masses); this translates into enhanced capture probabilities for equal masses. 



the gateways in the zero- v elocity surface at time t = (see 
Fig. were comparable lAstakhov fc Farrellvl l20o3) . Ap- 
proximately, this equates the fluxes of incoming/outgoing 
particles and facilitates a fair comparison of relative densi- 
ties of chaotic versus hyperbolic orbits. The actual energies 
used at t = (/ = 7r/2) were E — —1.725 for m r — 1 and 
E — —1.675 for m r = 0.05. Admittedly this is an ad hoc ap- 
proach but we find that the general appearance of the FLI 
plots is quite robust to the precise values of the initial Ja- 
cobi constant used provided that the gateways have similar 
sizes. 



4-3.2 Interpretation of FLI maps 

Figure |5] shows FLI maps for different mass ratios and ec- 
centricities. It is easy to distinguish between chaotic regions, 
directly scattering (hyperboli c) regions and - to the intru der 
- inaccessible KAM regions (lAstakhov fc Farrellvlliool . In 
the chaotic regions nearby orbits tend to be dense and have 
similarly large FLI values. This results in chaotic regions 
having the appearance of a shotgun-like pattern - unlike in 
a SOS, however, each point corresponds to the initial con- 
dition of a single orbit coloured according to the final com- 
puted value of its FLI. By contrast, in hyperbolic regions 
orbits rapidly leave the Hill sphere and enter heliocentric 
orbits. We chose 20 mutual orbital periods of the primaries 
(in the ERTBP limit) as a cut-off such that orbits which 
survived for less than this cut-off were not plotted. We ver- 
ified that our results are not sensitive to the precise choice 
of cut-off time. The large solid regions in Fig. |S| correspond 
to regular (KAM) regions. Elsewhere we have examined the 



effect of ellipticity on chaos-assisted capture and the differ- 
ences in FLI maps for chaotic, hyperb olic and regular zones 
was e xamined in detail in the ERTBP 1 Astakhov fc Farrellvl 
12004 . 

In summary, in FLI maps hyperbolic regions manifest 
themselves by their relative s parsity of points as compa red 
to chaotic and regular zones lAstakhov fc FarreUvlkoO^ . 

4-3.3 Effect of mass ratio and eccentricity on capture 
probabilities 

First we compare the situation for equal versus unequal 
masses in the case of circular orbits, i.e., Fig. [fji and 
HJa. Clearly the observable sea of chaos is much larger for 
frame (a) i.e., equal masses, than for frame (b) i.e., un- 
equal masses. However, in both cases it is also true that 
there are relatively few hyperbolic regions. For the case of 
unequal masses the regular KAM regions are very large 
but since KAM regions represe nt barriers in phase space 
iLichtenberg fc Liebermanll992l) these regions are inaccessi- 
ble to the intruder. Thus in both Fig. [(J, and |S|d an intruder 
is likely to become caught up in a chaotic zone thereby pro- 
ducing a resonance. According to our earlier arguments this 
will, on average, increase dwell-times and so tend to desta- 
bilize the binary. This suggests that circular (or very low 
eccentricity) binary orbits will be less likely to be captured. 

Next we examine the case of higher eccentricity, e = 0.7, 
shown in Fig. and |SJi. This case is actually the most 
relevant to the real situation since the osculating eccentric- 
ity of the chaotic binary is generally quite high. For equal 
masses, shown in frame (c), very large hyperbolic regions 
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are clearly visible whereas in frame (d) (unequal masses) 
the "non-regular" regions are primarily chaotic rather than 
hyperbolic. Incoming intruders will therefore encounter very 
different environments between the two cases. If the binary 
partners are of equal mass, intruders will most likely enter 
hyperbolic regions and so will exit the Hill sphere promptly. 
In the case of unequal masses the battleground is of a very 
different nature, resembling a field of sticky mud in which in- 
truders rapidly become bogged down. That is, they become 
temporarily trapped in trinary resonances which increases 
the probability that the binary will be disrupted rather than 
stabilized. This can be explained, in part, by noting that in 
the CRT BP with equal mass primaries (the Copen hagen 
problem, ISzebehelvl Il967t iBenet et alJll997t iNaeleJ l20oj) . 
the zero velocity surface is symmetric in x and so the gate- 
ways into and out of the capture region are equivalent. This 
facilitates quick transits through the capture zone, whereas 
in the asymmetric mass system one gateway is smaller than 
the other constituting a bottleneck which hinders passage. 

The FLI maps in Fig. |5] provide a qualitative explana- 
tion of the observed mass effect. The dwell-times of intruders 
within the Hill sphere are a sensitive function of whether in- 
truders enter hyperbolic or chaotic regions of phase space 
inside the Hill sphere. As system parameters change (total 
initial energy, ellipticity, etc.) the sizes and shapes of these 
regions alter in a highly nonlinear manner (see Fig. |3J;) as 
is characteristic of chaotic scattering. With increasing ellip- 
ticity the sizes of the chaotic versus the directly scattering 
regions decrease, especially in the case of equal mass binary 
partners. This translates into enhanced capture probabilities 
for equal masses. 

Overall, based on a comparison of the four cases in 
Fig. ©in the ERTBP, we conclude that intruder stabilization 
will be most efficient for (i) equal masses and (ii) moderate 
to high eccentricity orbits. It is somewhat ironic that, in this 
mechanism, the existence of chaos in the binary orbit (pro- 
duced by solar tides) is critical to nascent KBB formation, 
whereas if intruders enter zones of chaotic (as opposed to 
hyperbolic) scattering then the result is destabilization of 
the binary. 

However, the ERTBP model has its limitations. The 
assumption that the mutual chaotic binary orbit can be 
treated (from the point of view of the intruder) as an el- 
liptical Kepler orbit on short timescales is not valid if the 
actual binary orbit is itself extremely unstable on compara- 
ble timescales. Such three-body orbits tend to involve "near 
misses" of the binary partners and are less stable than are 
more "circular" three-body orbits (e.g., see Fig. |2K)- Thus 
for very high instantaneous eccentricities (orbital segments 
involving near misses) the binary partners (in the full four- 
body problem) are difficult to stabilize. This tends to re- 
duce the capture probability into very high eccentricity or- 
bits considerably. 

The net effect of all of these considerations is to favor 
moderately elliptical binary mutual orbits when the binary 
partners have approximately equal masses. This is precisely 
what is fo und in our simulations and in actual observations 
of KBBs <BurnJl2004lNoll et alJl2004ajbh . 



5 ALTERNATIVE FORMATION MODELS 

Of the four paths mentioned in the Introduction and sum- 
marized in Table extensive numerical simulations have 
only been reported for the exch ange mechanism, Path 4 
( Funato et al. 2004; Kolassa 2004). In this section we present 
a brief comparison of our multiple scattering model with al- 
ternative theories and also report some m odel simulations 
which we have performed for Paths 2 and 3 (iGoldreich et alJ 
2002). 



5.1 Exchange reactions 

In the exchange reaction model jFunato et al l2004h the first 
step is the formation of a binary through a two-body dis- 
sipative encounter. This can happen in two ways; (i) tidal 
disruption of an object during a close encounter followed 
by accretion of the resulting fragments to produce a bi- 
nary; and, (ii) a "giant" impact in which coagulating de- 
bris after the collision produces a satellite orbiting a larger 
object. Main-belt asteroid binar i es (and the Earth-Moon 
binary JCanup fc Asphaud 1200 ll : ICanud l2004h and (with 
some important d ifferences) possibly the Pluto-Charon bi- 
nary iCanurJ2005h 'l are thought to have formed in this man- 
ner. Generally, the result is a binary with an extreme mass 
ratio and a tight, circular orbit. In the exchange reaction 
mechanism, the mass ratio of such proto-binaries in the KB 
is increased through later gravitational scattering encoun- 
ters wit h a third object ori ginating (in the simulations) at 
infinitv. lFunato et al] ]2004f) modelled this by performing ex- 
tensive scattering simulations between strongly asymmetric- 
mass binaries already following bound, compact, circular 
orbits and large incoming intruder masses. They found a 
propensity for the smaller member of the initial binary to 
be replaced by the intruder with the final binary then fol- 
lowing a large, highly eccentric Keplerian orbit. 

While this mechanism is plausible and does produce 
binaries which are at least qualitatively similar to known 
KBB orbits a number of difficulties nevertheless remain. In 
the first place, this mechanism does not provide an ab initio 
explanation for why the the mass ratios of actual KBBs are 
of order unity. Instead, it simply provides a route by which 
a binary can increase its mass ratio. Because Funato et al. 
assumed that the primary member of the initial binary and 
the intruder had equal masses then, if exchange occurs, a 
binary of order unity mass ratio is guaranteed to result. 
There appears to be no co mpelling reason, at least based on 
the simulations reported jFunato et al]l2004h . to conclude 
that mass ratios of order unity should necessarily result from 
this process, only that exchange reactions can increase mass 
ratios. 

In addition, exchange reactions seem to produce orbital 
eccentricities ( and semimajor axes) larger t han those typi- 
cally observed <Nolil200St (Noll et aljl2004rj| . However, this 
might si mply be the re s ult of the particular choice of masses 
made bv lFunato et al] J2004) . Thus it would be interesting 
to extend these simulations to include a broader range of 
intruder masses so as to understand (a) the role of intruder 
mass in determining the probability of exchange and (b) the 
effect of intruder mass on post-exchange orbital properties 
of the binary. 
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5.2 Two-body collisions inside the Hill sphere 

The first mechanism propos ed for the production of trans- 
neptunian binaries is that of lWeidenschillinel l|2002l) . In this 
model two bodies accrete after a mutual collision within the 
Hill sphere of a third, larger, body. They then remain bound 
as a single object orbiting the larger body, thereby produc- 
ing a binary. The main objections to this model have to 
do with t he scarcity of la r ge obj ects in th e primordia l KB. 
However, iPetit fc Mousisl \20o4\ (see also ISternll2002l) dis- 
cu ss this issue at length and conclude that the mechanism 
of lWeidenschillind |.2002) might have operated, possibly in 
tandem with other mechanisms. In the absence of detailed 
simulations it remains an open question whether this mech- 
anism can provide quantitative agreement with the observed 
properties of KBBs. 

5.3 Capture through dynamical friction 

This m echanism was originally proposed bv lGoldreich et all 
(2002). Here we have shown that th e dynamical basis for 
this model is chaos-assiste d capture (lAstakhov et al.ll2003t 
lAstakhov fc Farrellvll2004) . No detailed simulations of KBB 
formation or t heir orbital and compos itional properties were 
presented by iGoldreich et alJ J2002l> . Because this mech- 
anism proceeds from similar assumptions to the one we 
propose, we performed a limited number of simulations in 
which multiple intruder scattering as an energy loss mecha- 
nism was replaced by dynamical f riction, modelled as simple 
velocity dependent dissipation jMurrav fc Dermotd Il999t 
lAstakhov et alJl2003l) . We find that, dynamical friction pro- 
vides an efficient route to binary capture and Keplerization, 
producing final orbital elements similar to those in Fig. [I] 
However, it is not sensitive to the mass ratio of the initial 
quasi-bound binary. This is because the scaled Hill equations 
0, even including dissipation, contain no parameters and so 
cannot depend on the masses of the binary partners. Thus, 
no mass effect is predicted. Of course, more sophisticated 
models of dynamical friction llC handrasekharlll 943l: iBinnevI 
Il977l: IGoldreich et al]l2004t 1 Just fc Penarrubiall2004l) might 
conceivably reveal such a preference. In addition, dynami- 
cal friction has been ruled out on the basis of estimates of 
planetesimal ma ss distributions alth ough this issue may be 
worth revisiting jFunato et al.ll20o"3) . 

5.4 The L 3 mechanism 

Agai n, this mechan ism was originally proposed by 
IGoldreich et al.l i2002T) and is the most similar to the one 
which we present here. A l thoug h no simulations of L 3 were 
made bv IGoldreich et ail <|2002T) a number of issues can be 
raised in regard to it; (a) in L 3 , as originally proposed, 
roughly equal mass binaries result directly from the assump- 
tion that a transient binary undergoes a scattering encounter 
with a third body of comparable mass to the binary part- 
ners (which are assumed to be themselves of similar mass). 
That is, the preference for order unity mass ratios is an 
inescapable result of the foundational assumptions of the 
model; (b) unless multiple L 3 type encounters are being as- 
sumed, then hardening and Keplerization of the large tran- 
sitory binary must occur in a single encounter which seems 
unlikely. To test this we have performed simulations exactly 
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Figure 7. Histogram showing the survival probability of equal 
mass chaotic binaries (ot\ = oti) after N enc L 3 encounters with 
large intruders (03 = 0.3, see equation JSJ). 1000 initial binaries 
were taken with initial conditions for three-body Hill orbits and 
intruders originated as described in Sec. |3] To generate statistics 
the simulation was repeated 50 times, each time choosing the 
phase of each binary randomly with respect to the intruder. 

as for Fig. [I] except using a large intruder. The results are 
summarized in Fig. [7] and indicate that multiple large body 
encounters tend to disrupt the binary catastrophically. Fur- 
ther, examination of the actual orbital elements produced in 
these simulations confirms that neither a single nor a small 
number of L 3 -type events is usually sufficient to produce 
Keplerian orbits. Instead, after a small number of such en- 
counters the orbits resemble those of the initial quasi-bound 
binary, e.g., that shown in Fig. [5] In these cases Keplerian 
orbital elements cannot meaningfully be defined. 

Of course a hybrid of our model and L 3 is also possible. 
Certainly our simulations do not rule out the possibility of 
L 3 -type collisions with additional hardening being produced 
through collisions with smaller bodies; however, L 3 collisions 
are not necessary for capture or to produce roughly equal 
mass binaries. 



6 CONCLUSIONS 

We have presented what are, to our knowledge, the first 
simulations of binary-intruder scattering in which the bi- 
nary orbit is itself chaotic. Previous simulations of binary 
star intruder scattering have focussed on the important but 
simpler case in which the initial binary already follows a 
bound Keplerian ellipse. Here the binary orbits used are 
long living chaotic orbits (in the Hill approximation) which 
cling for long periods to regular KAM structures in phase 
space. These transitory orbits are possible only if the Sun- 
binary Hill sphere is taken into account. The importance 
of the Hill sphere in tem porarily capturing Kui per-belt bi- 
naries was first noted bv lGoldreich et alJ (l2002f) on empiri- 
cal grounds. The dynamical explanatio n for the results pro- 
posed here is chaos-assisted capture lA stakhov et al-lEoOSl : 
lAstakhov fc Farrellvll2004l: iHolman et al.ll2004li . 

Once a proto-binary has formed within the Hill sphere 
then, if it is to survive, some stabilization mechanism must 
operate. Our simulations demonstrated that multiple low- 
mass intruder scattering can not only stabilize transitory 
binaries but also provides an efficient route to binary hard- 
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ening and Keplerization of the orbit. By adjusting a sin- 
gle parameter - the number of intruder-binary encounters - 
we were able to obtain excellent agreement with the orbital 
properties of known KBBs. In particular, our simulations 
predict that KBBs will have moderate eccentri cities . This 
is in good agreem ent with recent observations l)Nollll2003t 
iNoll et al.ll2004allrl and contrasts with the exchange mech- 
anism of KB B formation which produces exclu sively large 
eccentricities (F unato et al 1 l2004tlKolassall200^ . 

Our simulations also lead to a striking preference for 
the production of binaries having roughly equal masses. 
Strongly asymmetric mass ratios - typical of main belt as- 
teroid binaries - are selected against, and, again this appears 
to be consistent with observations. The model we propose 
is not set up to produce only order unity mass ratio bina- 
ries; nevertheless, a preference for mass ratios of order unity 
emerges. 

An important difference between our model and L 3 con- 
cerns how the binary hardens. Single encounters with large 
intruders are unlikely to produce the relatively compact (as 
compared to Rh) Keplerian binary orbits that are actually 
observed. Moreover, we find that multiple encounters with 
large intruders tend, on average, to disrupt the binary. In- 
stead, the model we propose involves a sequence of low-mass 
intruder events rather than a single scattering encounter 
with a massive object. It turns out that this provides an 
efficient mechanism for binary hardening and Keplerization. 
However, because we are unable to estimate the probabili- 
ties for the initial rate of production of quasi-bound binaries 
of varying mass ratios in the Hill sphere, it is unclear how 
strongly this mass effect will persist. Nevertheless, our simu- 
lations demonstrated that equal mass binaries have a much 
higher survival probability after multiple encounters than 
do asymmetrical mass ratio objects. Thus, we are confident 
that the mass effect we have found will not be washed out 
by the statistics of transitory binary formation. 

Throughout we have assumed that the binary centre- 
of-mass follows a circular (heliocentric) orbit. However, in 
principle, we also need to consider the effect of introduc- 
ing eccentricity into the heliocentric binary orbits. This is 
because some observed KBBs have significantly eccentric he- 
liocentric orbits with perihelia close to, or even inside, Nep- 
tune's orbit (e.g., Pluto-Charon, (47171) 1999 TC 36 and 
(26308) 1998 SMws)- Typical KB B heliocent ric elliptici- 
ties lie in the range ~ 0.03 - 0.37 jNol]||2003fi . In our re- 
cent study of captur e and escape in the elliptic restricted 
three-body problem (lAstakhov fc Farrellvl 120041) we found 
that the introduction of moderate heliocentric ellipticity re- 
duces, but does not entirely eliminate, the efficacy of the 
chaos-assisted capture mechanism. While further study of 
this issue is clearly warranted, we expect that the results 
presented here will similarly carry over to the elliptic helio- 
centric case. 

Of course, it is also possible that observed heliocen- 
tric TNO binary eccentricities are not primordial. Actually, 
this is suggested by the observation that massive bodies in 
eccentric orbits will tend to have large approach velocities 
apart from coincidences such as, e.g., when perihelia are 
aligned. Because such cases are rare, relative velocities will 
tend to be large and this will serve to reduce the probabil- 
ity of capture into chaotic orbits. This conclusion is consis- 
tent with results for capture in the elliptic restricted three- 



body problem (lAstakhov fc Farrellvll2rjrj4T) . These consider- 
ations suggest that, as with Paths 1-3, binaries formed 
when the disk was massive and dynamically cold, i.e., prior 
to t he dynamical excitation and depletion of the Kuiper- 
belt jDuncan fc Levisonll997l : lGomesl2003l iMorbidelli et alJ 
l2004h . 

Finally, we note that Figs. 13151 demonstrate that, while 
significant, the selection for order unity mass ratio objects 
is not to the complete exclusion of smaller mass ratios. For 
example, the differences in dwell-time distributions are not 
so great as to eliminate the capture of small mass-ratio bi- 
naries entirely. Rather, Figs. [3] and 2] lead us to predict that 
a sizeable population of asymmetric mass KBBs is prob- 
ably awaiting detection. Heralds of this group of objects 
might already exist in the low mass ratio binary objects 
(4717 1) 1999 TC 36 JTruiillo fc Brownll2002l: lAltenhoff et alJ 
| 2004h and (26308) ^998 SM^ jBrown fc Truiillol 120021 : 
ISternll2002l: lNolll2003l : IPetit fc Mousisll2004h . 
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